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Bragg scattering of light in a strongly interacting trapped Fermi gas of atoms 
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Physical Research Laboratory. Navrangpura, Ahmedabad 380 009, India 
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We study Bragg scattering of laser light by trapped Fermi atoms having two hyperfine spin 
components in the unitarity-limited strongly interacting regime at zero temperature. We calculate 
the dynamic structure function of the superfluid trapped Fermi gas in the unitarity limit. Model 
calculation using local density approximation shows that, the superfluid pairing gap in the unitarity 
limit is detectable from the measurements of dynamic structure function by Bragg spectroscopy, 
while in the weak-coupling BCS limit, the gap eludes such spectroscopic detection. 

PACS numbers: 03.75.Fi,74.20.-z,32.80.Lg 



I. INTRODUCTION 

Achieving Bardeen-Cooper-Schrieffer (BCS) type 
phase transition in trapped neutral fermionic atoms is 
the primary goal in the current experimental research 
with ultracold atomic Fermi gases. Recently, two groups 
0, 0] have measured collective oscillations in strongly 
interacting trapped Fermi atoms. The results of their 
measurements suggest the occurrence of fermionic super- 
fluidity in the atoms according to a theoretical prediction 
0. Condensation of "fermionic atom pairs" has been 
reported earlier Q. The superfluid pairing is believed 
to occur near the crossover yj |(| between the predicted 
BCS state of fermionic atoms and the Bose-Einstein con- 
densation of bosonic dimer molecules formed from Fermi 
atoms due to a magnetic field Feshbach resonance. Sev- 
eral groups have experimentally produced Bose-Einstein 
condensates (BEC) of molecules formed from degen- 
erate Fermi atoms in the vicinity of the Feshbach res- 
onance. The atoms become strongly interacting at and 
near the Feshbach resonance. Several other recent experi- 
ments 0, 0, 0, ^3 with Fermi atoms have demonstrated 
the strong-coupling behavior of the atoms. Therefore, 
many-body effects become important in describing the 
physics near the predicted BCS-BEC crossover in Fermi 
atoms. Recently, a number of theoretical investigations 
|12| have revealed many intriguing aspects of many-body 
effects in the crossover regime. Despite the recent experi- 
mental indications of the occurrence of Cooper-like pair- 
ing near the crossover, the energy gap for such pairing 
has not been so far measured due to the lack of efficient 
detection methods. 

There are several methods available for measuring pair- 
ing gap in electronic superconductors, but any of them 
is hardly applicable to the trapped neutral Fermi atoms. 
Recently, several theoretical proposals [lH. IbH ] have been 
made for detecting atomic Cooper pairs by laser light. 
There has been a suggestion [l3| to use resonant light 
for exciting one of the Cooper-paired atoms into an ex- 
cited electronic level, and thus to make an interface be- 
tween normal and superfluid state in analogy with the 
well known tunnelling experiment in superconductors. 
Our purpose here is to examine Bragg scattering of lasers 
as a method for measuring the paring gap of a super- 



fluid trapped Fermi gas of atoms. Bragg spectroscopy 
has been already used for measuring structure function 
of an atomic BEC 15] . Unlike in a BEC, the analysis 
of Bragg scattering in a superfluid trapped atomic Fermi 
gas is complicated due mainly to the Pauli blocking and 
the paring gap. 

In this paper, we study Bragg scattering of light in a 
two-component trapped atomic Fermi gas at zero tem- 
perature. Bragg scattering in a Fermi gas of atoms has 
an analogy with Raman scattering in electronic supercon- 
ductors. We therefore, develop our theoretical treatment 
followin g th e theory of Raman scattering in supercon- 
ductors The relevant physical quantity is the dy- 
namic structure function which is the Fourier transform 
of density-density correlation function in time 17]. It is 
a measure of the spectrum of density fluctuation and pro- 
portional to the differential scattering cross section per 
unit energy. In the case of superconductors, it is well 
known that unless the energy transferred from a probe 
to the superconductor exceeds 2A (A being the pairing 
gap), the dynamic structure function is zero. At 2A, it 
shows a sharp discontinuity. In the case of superfluid 
trapped Fermi gas, A has a spatial distribution vary- 
ing from a maximum at the trap center to a vanishingly 
small value at the edge of the trap. We find that the dy- 
namic structure function of a superfluid trapped Fermi 
gas in the unitarity limit has a prominent shift compared 
to that of a normal or superfluid trapped Fermi gas in 
the BCS limit. It also exhibits a discontinuity at 2A(0) 
where A(0) refers to the energy gap at the trap cen- 
ter. As the energy transfer decreases below 2A(0), the 
dynamic structure function falls off with reducing slope, 
and as the energy transfer goes to zero, the slope van- 
ishes. In contrast, in the case of a normal trapped gas 
or a superfluid trapped gas in the BCS limit, the dy- 
namic structure function increases almost linearly with 
the energy transfer in the low energy regime. 

The paper is organized as follows. In the following sec- 
tion, we discuss in some detail the physical scenario for 
Bragg scattering in two-component trapped Fermi atoms 
in the unitarity-limited strongly interacting regime. In 
Sec. Ill, we discuss the theoretical treatment for calculat- 
ing dynamic structure function of a uniform superfluid 
Fermi gas. We then generalize this treatment for a su- 
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perfluid trapped Fermi gas in Sec. IV. The results are dis- 
cussed in Sec.V. We conclude in Sec. VI. 



II. TRAPPED FERMI GAS IN THE 
UNITARITY REGIME 

We consider a harmonic trap with potential Vho(r, z) = 
(l/2)m(u>±r 2 _ + u) z z 2 ), where lu± and u> z denote the ra- 
dial and axial frequency, respectively, of the trap. The 
harmonic oscillator is characterized by the radial (axial) 
length scale au z ) — •JhJ (nuo±M. One can define a 
geometric mean frequency of the harmonic oscillator by 
uj} lo = (uj 2 _uj z ) 1 /' i and a corresponding geometric mean 
length scale of the oscillator by ay lo — ^h/(muJho)- In 
our treatment, we resort to Thomas- Fermi local density 
approximation (LDA) |18| which is particularly applica- 
ble when the local Fermi energy is larger than the average 
level spacing of the trap and the coherence length of the 
fcrmion-pair is shorter than the average trap size. Under 
this approximation, the state of the system is governed 
by 



e F (r) + V ho (r) + U(r) = /z 



(1) 



where ep(r) = h 2 kp(v) 2 / (2m) is the local Fermi energy, 
k F (r) denotes the local Fermi momentum which is related 
to the local number density by n(r) — kp(r)' i / (Qtt 2 ). 
Here U represents the mean-field interaction energy and 
fi is the chemical potential fixed by normalization condi- 
tion. At low energy, the mean-field interaction energy 
depends on the two-body s-wave scattering amplitude 
fo{k) = —a s /(l + ia s k), where a s represents energy- 
independent s-wave scattering length and k denotes the 
relative wave number of two colliding particles. In the 
dilute gas limit (|a s |fc « 1), U becomes proportional to 

a s in the form U(r) = iv 2n ^ s n ( r )- I n the unitarity limit 
\a s \k — > oo, the scattering amplitude /o ~ i/k and hence 
U becomes independent of a s . It then follows from a sim- 
ple dimensional analysis that in this limit, U should be 
proportional to the Fermi energy: U(r) = (3ep(r) where 
(3 is the constant of proportionality. Under LDA, the 
density profile of a trapped Fermi gas is given by 



n(r) = n(0) 



R 2 



3/2 



where 



n(0) = 



1 



2m M \ 3/2 
f>V/Y ; \1 + (3J 



(2) 



(3) 



is the density of the atoms at the trap center. Here 
R\i z \ — 2/V( mw l(z)) being the radial(axial) Thomas- 
Fermi radius. The normalization condition on Eq. @ 
gives an expression for 



where N a is the total number of atoms in the hyperfine 
spin a. The Fermi momentum kp = [37r 2 n(0)] 1 / 3 is then 
given by 



k, 



where 



, _ (48V CT ) 1 /6 



aho 



(5) 



(6) 



H = (l+l3) 1/2 (6N a ) 1 / 3 f l io 



(4) 



is the Fermi momentum of the noninteracting trapped 
gas. For an attractive interaction (— 1 < j3 < ) in the 
unitarity limit, Fermi momentum would be larger than 
that of noninteracting gas by a factor of (1 + /3) -1 / 4 . 
Accordingly, the Fermi energy would be larger than the 
Fermi energy e F — hw F of noninteracting gas by a factor 

of {i + py 1 / 2 . 

To illustrate the main idea, we specifically consider 
trapped 6 Li Fermi atoms in their two lowest hyperfine 
spin states | g\) =| 2Si/ 2 ,-F = 1/2, m F = 1/2) and 
| g-i) =| 2Si/ 2 ,F = 1/2, nip — —1/2). For s-wave pairing 
to occur, the atom number difference SN of the two com- 
ponents should be restricted by ^ < T c /ep where T c is 
the critical temperature for superfluid transition and tp 
is the Fermi energy at the trap center. Unequal densities 
of the two components result in interior gap superfluidity 
[Ifi Efil ] which naturally arises in QCD matter. For sim- 
plicity, we consider the case N 1 / 2 — N_x/2 which is the 
optimum condition for s-wave Cooper pairing. An ap- 
plied magnetic field tuned near the Feshbach resonance 
(~ 822 Gauss) results in strong inter-component s-wave 
interaction. At such high magnetic fields, the splitting 
between the two ground hyperfine states is ~ 75 MHz 
|2l| . while the corresponding splitting between the ex- 
cited states | ei) =| 2P 3/2 ,.F = 3/2, m F = -1/2) and 
I e 2 ) =| 2P 3/2 , F = 3/2, m F = -3/2) is - 994 MHz J^. 
Taking advantage of these Zeeman splits, it is possible 
to scatter atoms selectively of one hyperfine spin compo- 
nent by making use of a pair of circularly polarized laser 
beams-both having either <j+ or cr_ polarization. 

In Bragg scattering, two laser beams (Bragg pulses) 
with a small frequency difference are applied to the 
trapped gas. The magnitude of this momentum trans- 
fer is q ~ 2kL sin(6>/2), where 9 is the angle between 
the two beams and kp is the momentum of a laser pho- 
ton. Let both the laser beams be ct_ polarized and tuned 
near the transition | 172) — H £2)- Then the transition 
between the states | g\) and | e 2 ) would be forbidden 
while the transition | gi) — >| ei) will be suppressed due 
to the large detuning ~ 900 MHz. This leads to a sit- 
uation where the Bragg-scattered atoms remain in the 
same initial internal state | e 2 ). Similarly, atoms in 
state I g\) only would undergo Bragg scattering when 
two 0+ polarized lasers are tuned near the transition 
I g x ) ->| 2P 3/2 ,F = 3/2, m F = 3/2). Thus, we infer 
that in the presence of a high magnetic field, it is possi- 
ble to scatter atoms selectively of either spin components 
only by using circularly polarized Bragg lasers. 
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For simplicity, let us first consider Bragg scattering in 
a uniform Fermi gas. Later, we will analyze the nonuni- 
form case. We assume that both the laser beams are ct_ 
polarized and tuned near the transition | c/2) — H The 
detuning of both the lasers from the atomic resonance 
should be much greater than the single- and two-photon 
line widths in order to avoid heating of the system. Un- 
der such conditions, for a uniform gas, the laser-atom 
interaction Hamiltonian can be expressed as 



along with the equation 



H, 



k 



c{(k + q) Ci (k) + H.c. 



(7) 



where Cj (k) represents the atomic field operator with mo- 
mentum k, the subscript J. (t) refers to the state | (72) 
(| <?i)) and is the two-photon Rabi frequency. Here 
5 = lo\ — LU2 is the detuning between the two lasers. One 
can identify the operator /oj(q) = 5TJ k c[(k + q)c|(k) as 
the Fourier transform of the density operator pi(r) — 
^j(r)^(r) where &j(r) represents the field operator in 
the real space. The spectrum of the scattered atoms 
would be proportional to the rate of transition proba- 
bility which, according to Fermi's Golden rule, is given 
by 

K = 2^ 2 ^|(/|p t i (q) \0)\ 2 S(HS~e f + e Q ) 
f 

= 27rM} 2 S(q,w) (8) 

where |0) represents the many-body ground state and 
sum runs over all the final states |/) which can be coupled 
to the ground state by the operator p{. 



III. DYNAMIC STRUCTURE FUNCTION OF A 
UNIFORM SUPERFLUID 



Using Bogliubov transformations 

C|(k) = Ufe7i(k) + Vfc7j(-k) 



c T (k) = Mfe7 T (k) - Ufe7;(- k ) 



(9) 



(10) 



one can reexpress the interaction Hamiltonian in terms 
of quasiparticle operators 7^(k) and 7 CT (k). Here the nor- 
malization condition is lu^-l 2 + \v k \ 2 — 1 with 



Kl 2 = (i/2)(i 



and 



Ek = je k +Ai 



(11) 



(12) 



where = h 2 k 2 / (2m) — /i and A k is the pairing gap. The 
chemical potential fi and the gap A k can be obtained by 
solving the regularized gap equation 



\ " 



V^\2e k 2E k ) 



(13) 



— e -iv 

k 



E k 



(14) 



of the density of single component. Here V is the vol- 
ume of the system. The analytical solutions of these two 
coupled equations have been already obtained in the Ref 
23] for the entire range of the parameter a s k F starting 
from weak interaction (a s k F — > ±0) to the unitarity limit 
(a s k F — > ±00). In the unitarity limit, the mean-field in- 
teraction becomes independent of the scattering length 
and the Fermi system exhibits universal behavior. In 
this limit, the only energy scale available to the system 
is the Fermi energy. The mean-field energy and the gap 
then become proportional to the Fermi energy. Solutions 
of the above two equations in the unitarity limit yield 
A ~ 1.16/i, fJ, = (1 + P)e F , where = -0.41 [H H3 
is a constant. Calculations using quantum monte carlo 
simulation [25| yields (3 — —0.56 and A ~ 0.49ei?, while 
calculations using Galitskii and lowest order constraint 
variation (LOCV) approximations j2|| give f3 = —0.67 
and P = —0.43, respectively. The constant (3 was experi- 
mentally introduced and measured to be -0.26 in Ref 
while it was experimentally found to be between -0.4 to 
-0.3 in Ref 01 ■ A recent experiment [27| has obtained 
(3 = -0.68. 

One can notice that the BCS ground state |0), which 
is annihilated by the quasiparticle operators, will have 
nonzero matrix element of p q only when the final state 
I/) is a state of two quasiparticles with mutually opposite 
spins and their momenta differing by q. The dynamic 
structure function thus takes the form 

S(q,S) = J2\ u k' v k\ 2 6(hS-E k , -E k ) 



V 



4(2tt) 3 
x 6(hS - E k , - E k ) 



ffl ^ Ek ' + £k')(E k - £fc) 



E k > E k 



(15) 



where k' = |k + q| is the wave vector of a scattered atom. 
Note that the usual BCS coherence factor m(k, k') = 
u k >v k + u k v k i 0,13 which appears in the description of 
electronic superconductors has changed. This is due to 
the fact that the polarization-selective dipole transitions 
in Fermi atoms in the presence of strong magnetic field as 
discussed earlier lead to the transfer of momentum and 
energy to either partner (of hyperfine spin f or I) of a 
cooper pair, while the other partner remains almost im- 
mune to the momentum and energy transfer. However, 
since a particle's state is a superposition of two quasipar- 
ticle states of both the spins, both the spin states will be 
affected in the quasiparticle framework . The presence 
of the (5-function in the integrand reveals that S*(q, S) 
would be nonzero only when HS > 2A. The integration 
should be carried out subject to the restrictions k < kp 
and k' > k F resulting from Pauli blocking. Following 
the method of Ref 0, we evaluate the integral. After a 
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lengthy algebra (the method of calculation is described 
in the appendix), one can express the Eq. Q15[) in the 
form 



SfaS) 



vp A 2 



2 p q VFhS 
x Fdz (1+JN)2 

Jo (i-i^ 2 ) 3 / 2 (i -^ 2 ) 1 /2' 



(16) 



where i/p is the density of states at the Fermi surface. 
P q =frq,j = l- 4A 2 /(M) 2 and 



z Q = Min 



1, 



PqVF 



(H 2 8 2 - 4A 2 )V2 



(17) 



If 2A < M < (p q v F ) 2 + 4A 2 ) 1 / 2 , then z = 1 and the 
result is 



at x\ Vp S 
SfaS) = — 

8 qvF 



E(j) + J -{8 + M 2 F 1 (3/2,3/2;2J)} 



(18) 



where M = —j). Here E(j) represents the complete 
elliptic integral and 2-Fi(a, b; c, d) is the hypergeometric 
function. In the limit A — > 0, S(S,q) — vpS / '(2g«j?) 
which is exactly the dynamic structure function of normal 
fluid H3. For HS > (p q v F ) 2 + 4A 2 ) 1 / 2 , SfaS) can be 
expressed in terms of the Elliptic integrals of first and 
second kind. In this case, as A — * 0, Sfa S) — ► and 
becomes independent of q. Since this case is not suitable 
for obtaining information about the gap, we henceforth 
focus our attention only on the former case, that is, 2A < 
HS< (p g w F ) 2 +4A 2 ) 1 /2. 



IV. DYNAMIC STRUCTURE FUNCTION OF A 
TRAPPED SUPERFLUID 

Let us now turn our attention to the dynamic structure 
function of a superfluid trapped Fermi gas. This can 
explicitly be written as 

SfaS) = I / d 3 x<(xK„(x)exp(iq.x) | 2 



x 5(S-E n -E m ) 



(19) 



In the case of a uniform gas, we have Uk(x) = 
Uk exp(ik.x) and u_k(r) = v k exp(— ik.x), and hence the 
dynamic structure function defined by Eq. (|19fl takes the 
form of Eq. i|15|) ■ In the LDA, it can be expressed as 

Slda{<1,5) = J d 3 x J d 3 kn fe ,(x)[l - n fc (x)] 

x 5(5-E k ,(x)-E k (x)) (20) 

where n&(x) is the local momentum distribution defined 
by 

n fc (x) = |t. fc (x)| 2 =hl- ffc(x)/£? fc (x)]. (21) 



Here ^fe(x) = h k : /(2m) — /x(x), /i(x) = (J, — V?j (x) and 
i?fc(x) = -\/^fc(x) 2 + A(x). The momentum distribution 
of trapped Fermi atoms has been already calculated in 
Ref |24| . After performing the integration over k as in 
the preceding section, one can obtain 



Slda (q, 5 ) = \J d 3 xn(x) 
x E( Jx ) 



A(x) 



p q v F (x)tid~ 



f {8 + M 2 ^(||2, Jx )} 



(22) 



where u(x) = (3/2)n(x)/eir(x) is the local density of 
states per unit volume. Here j(x) = 1 — 4A(x) 2 /(M) 2 . 
The integration over the volume must be carried out 
subject to the boundary condition 2A(x) < HS < 
[{p q v F {x)) 2 + (2A(x)) 2 ] 1 / 2 . This means that, / d 3 x = 



rdr Io 



R z y/l-r2/R ± 



dz where r m i n = if 2A(r 



0, z = 0) < HS, otherwise it is the solution of the equation 
2A(r, z = 0) = HS. Here r max < R± is the solution of the 
equation [(p q v F (r, z = 0)) 2 + (2A(r, z = 0)) 2 ] 1 / 2 = M. 

In the BCS limit (fc F a s — > _ ), the gap is exponentially 
small and can be expressed by the well known formula 



8e F , 7r 
^bcs - — exp(- )■ 



(23) 



In calculating S(S, q) of trapped atoms in the BCS limit, 
one can follow the same method of calculations as de- 
scribed above. However, the parameter fejr | a s | in the ex- 
ponent of Eq. ll'-il) may be approximated by an average 
value. 



V. RESULTS AND DISCUSSIONS 

Figure 1 shows the results of our calculations. Plot- 
ted is the dynamic structure function S (8, q) of super- 
fluid trapped Fermi atoms as a function of energy trans- 
fer for different values of momentum transfer q. S(S,q) 
has been scaled by v% = (3/2)y TF (fc° ) 3 /(67r 2 e^), where 
Vtf — t^R\R z - In the unitarity limit, the behavior of 
S (<5, q) is quite different from that of superfluid in the 
weak-coupling BCS limit. This can be attributed to the 
occurrence of large gap in the unitarity limit. In plotting 
dashed-dotted curve corresponding to the BCS limit, we 
have taken A; F |a s | = 0.3, i.e., Abcs — 0.006e F . For this 
small value of the gap, S(8, q) reduces to almost that of 
normal fluid. We have checked this by calculating S(S, q) 
for a trapped normal fluid; and S(S, q) for normal fluid 
is almost indistinguishable from that of BCS superfluid. 
This can also be checked, as discussed in Sec. Ill, by tak- 
ing the limit A — > for which S(S, q) reduces to that of 
normal fluid. 

Particularly distinguishing feature of S(S, q) in the 
unitarity-limited strongly interacting regime as compared 
to normal or BCS regime is the shift of the peak. This can 
be noticed by comparing the dashed-dotted (BCS) and 
solid (unitarity) curves which are plotted for the same 
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FIG. 1: Dimensionless dynamic structure function S(S, q)/u F 
of a trapped superfiuid Fermi gas is plotted as a function of di- 
mensionless energy transfer S/lj f for the momentum transfer 
q/kF ~ 1.5 (solid), q/kF ~ 1 (dashed), q/kF = 0.5 (dotted) in 
the unitarity limit with f3 = -0.56 and A = 0.49e F = 0.74e F . 
For a comparison, S(S, c\)/v% in the BCS limit (dashed- 
dotted) with |a s |fci? = 0.3 is plotted for q/kF = 1.5. The 
shift for the unitarity limit (solid) compared to the BCS one 
(dashed-dotted) is 0.8A 

momentum transfer (q/k F = 1.5). This shift is propor- 
tional to the gap. For the parameters chosen in Fig.l, 
the shift is ~ 0.8A = 0.8 x 0A9e F = 0.8 x 0.49 x (1 + 
/3) _1 ^ 2 e^ — 0.6e F . It should be mentioned here that the 
dynamic structure function of a superfiuid trapped Fermi 
gas in the BCS limit has been evaluated earlier using a 
different method of calculation [29j which also shows no 
conspicuous shift except the appearance of an asymmet- 
ric peak. 

Another interesting feature is the discontinuity in 
S(S,q) at energy HS — 2A(0), where A(0) refers to the 
gap at the trap center. For a uniform Fermi superfiuid, 
S (S, q) remains zero until energy transfer exceeds 2 A at 
which it rises sharply with the increasing energy transfer. 
For a superfiuid trapped Fermi gas, owing to the spatial 
distribution of the gap, the dynamic structure function 
has a structure below 2A(0). As the energy transfer goes 
to zero, the gradient of S(S, q) vanishes. In the low en- 
ergy regime (HS < 2A(0)), S(S, q) varies with energy non- 
linear ly. When the energy transfer approaches to 2A(0), 
the gradient changes abruptly implying the discontinuity. 
This behavior can be explained by considering the bound- 
ary conditions 2A(x) < HS < [(p 9 i; F (x)) 2 + (2A(x)) 2 ] 1/2 - 
The lower bound on HS implies that, when HS is less than 
2A(0), the atoms at the central region of the trap can 
not respond to the Bragg pulses, only those atoms in 
the peripheral region can be responsive. As HS increases 
from the origin, increasing number of atoms from the in- 
ner region of the trap start responding. As HS — ► 2A(0), 
the atoms at the trap center can respond. On the other 
hand, due to upper bound on HS, as HS increases, in- 



creasing number of atoms from the peripheral region will 
cease to respond. Thus, as HS exceeds 2A(0), for large 
momentum transfer, S(S, q) will rise to a maximum and 
then eventually vanish as HS reaches its upper bound 
{(p q v F (0)) 2 + (2A(0)) 2 ] 1/2 - 

In the recent experiments 0, 0, |n| with two- 
component 6 Li atoms, the typical value of the Fermi ve- 
locity is vf = Hkp/m ~ 15 cm/second. The wavelength 
for transition 2Si/2 — > 2P 3 / 2 in Li atoms is A ~ 670.8 nm. 
With counter-propagating Bragg pulses tuned near this 
transition, the momentum transfer would be q ~ 2/c^x, 
where k F = 2ir/\. This momentum transfer raises the ve- 
locity of the scattered atoms by 2 x Hki,/m ~ 20 cm/sec 
which exceeds vp. Therefore, the scattered atoms should 
be distinguishable in time of flight images and hence 
dynamic structure function should be measurable. The 
polarization-selective Bragg spectroscopy as discussed in 
Sec. II may lead to betterprecision in time-of-flight spin- 
selective measurements 0^0 of the scattered atoms. Be- 
cause, the scattered atoms at low temperature would be 
almost noninteracting. 

Although, we have not studied dynamic structure func- 
tion in the BEC limit (a s kp — > + ) in which Fermi atoms 
form bosonic molecules which then condense to form a 
molecular BEC, the basic features of the dynamic struc- 
ture function in this limit can be understood intuitively. 
Unless the energy transfer exceeds the molecular binding 
energy e m = ft 2 /(ma 2 ), the molecules would not be dis- 
sociated into the atoms. For energy transfer greater than 
the binding energy, the dissociated atoms would have 
normal Fermi distribution, because they are repulsively 
interacting [a s is positive). Hence, dynamic structure 
function would be similar to that of normal fluid except 
the small shift e m which may not be even discernible. 



VI. CONCLUSION 

In conclusion, we have studied the Bragg scattering of 
light in a superfiuid trapped Fermi gas in the unitarity 
limit. Our results suggest that it is possible to detect 
the pairing gap in this limit by large-angle (i.e., large 
q) Bragg scattering. At small momentum transfer, the 
scattered atoms may not be distinguishable in the time- 
of-flight images. Hence, it may be difficult to observe 
experimentally all the features of the dynamic structure 
function in the low energy regime. In the weak-coupling 
BCS limit, the gap is so small that it would elude de- 
tection. The possibility of detecting pairing gap in the 
unitarity limit has been earlier discussed qualitatively in 
Ref [2(j. In this paper, we have provided quantitative 
justification of this possibility. More exact calculations 
should involve solving Bogoliubov-de-Gennes equations 
and the gap equation at finite temperature for the entire 
range of interaction strength. 



6 



Acknowledgments 

The author is grateful to G. S. Agarwal and S. Dutta 
Gupta for encouragement and discussions, and to M. 
Randeria for a helpful suggestion. He is thankful to P. K. 
Panigrahi, H. Mishra, A. Mishra and A. Dilip K. Singh 
for many stimulating discussions. 



APPENDIX A 

We here outline the method of calculation of the inte- 
gral in Eq. 1)15(1. The dominant contribution to the inte- 
gration comes from the fc-values near the Fermi surface. 
Therefore, restricting the integration near £ = efe — /i ~ 0, 
the Eq. H15fl can be reexpressed as 



S(q,6) 



V 
4(2tt) 3 
V 



J d 3 k J d£6(£ h -0F(q,Z) 
4(2tt)3 /^fo) SZ WW ( A1 ) 



where F(k, £) 
and E' = (£' 2 



(E'+j')(E-(,) 



— L with E = (e + A 2 ) 1/2 

A 2 ) 1 / 2 . It is of advantage to change 
the variable of integration into E by using the relation 
d£ = EdEj (E 2 - A 2 ) 1 / 2 . Considering E' as a function 
of E and defining the function F(E) = E + E'(E), one 
can use the identity 



5(u-E- E') = 



6(E - Eq) 
dF/dE \e=e 



where Eq is the solution of the equation F(E) = 0. After 
a lengthy algebra as in Ref 0] , one finally obtains 



V 



(27T) 



d 3 k5(£ k ) 



A 2 (M+| Pg .v fc | 



X> 3 / 2 X> 



■1/2 



(A3) 



where p q = hq, v k = Kk/m, V = h 5 — (p g .Vfe) and 
V = V- 4A 2 . Writing d 3 k = {2nm 2 /h 3 )de k v k sin 9 dO, 
where 9 being the angle between k and q, the integration 
over efe results in 



5(q,<S) = ^ f A 2 Jdx 



(H8 +p q v F \x\f 



(A4) 



where x — cosf? and Vp = (3N a /2ep) is the density of 
states of the single spin a. Here subscript F implies 
that the functions V and V are evaluated at the Fermi 
surface, that is, at v k — vp. Changing the variable of 
integration into 



PqVpX 



(fi 2 5 2 -4A 2 ) 1 /2' 



(A5) 



(A2) one can express the Eq. I|A4|I in the form of Eq. ijTfSJl 
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